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1.  Introduction 

This  paper  investigates  the  estimation  jf  the  parameters  (both  loca¬ 
tion  ar\d  scale)  of  the  logistic  distribution  using  sample  quantiles  and 
order  statistics.  Three  kinds  of  estimators  have  been  considered;  (l)  Best 
linear  unbiased  estimators  based  on  sample  quantiles;  (2)  Unbiased  linear 
asymptotically  best  estimators  of  Blom;  (3)  Asymptotically  best,  asymptotically 
unbiased  linear  estimators  of  Jung.  All  these  methods  of  estimation  are  asymp¬ 
totically  efficient  and  one  of  the  purposes  of  this  investigation  is  to  deter¬ 
mine  how  good  they  are  when  compared  to  the  best  linear  unbiased  estimators 
in  terms  of  their  relative  efficiency. 

Let  (x^,  Xg,...,  xQ)  be  a  sample  from  a  logistic  distribution  with 
p.d.f.  of  x  given  by 
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where 


6.  •»  n/  \/3~ 


The  c.d.f.  F(x)  is  then  defined  by 


(1.2) 


F(x) 


1 

!  +  e-a(x-H*Vo 


Section  2  contains  a  discussion  of  the  estimators  based  on  sample  quan¬ 
tiles.  When  a  is  known  the  optimum  symmetric  spacing  of  the  quantiles 
used  in  the  estimation  of  y  has  been  obtained,  for  ary  number  k  of 
quantiles.  When  |i  is  known,  the  optimum  spacing  of  the  quantiles  for 
estimating  CT  has  been  derived  for  k  =  3, 

Section  3  contains  a  brief  discussion  of  the  approximations  to  the  best 
linear  unbiased  estimators  suggested  by  Elcm  (1957)  and  Jung  (1956),  and  com¬ 
pares  these  two  sets  of  estimators  in  terms  of  their  relative  efficiency  with 
respect  to  the  best  linear  unbiased  estimators. 

2.  Quantile  Estimators 

The  quantile  estimators  are  based  on  a  fixed  number  of  sample  quantiles 
when  the  total  sample  size  is  very  large.  Such  a  method  of  estimation  would 
thus  be  useful  when  the  experimenter  has  a  very  large  sample,  but  would  like 
to  estimate  the  parameters  with  a  few  selected  observations  which  he  has  the 
freedom  to  choose.  Such  a  situation  could  arise  in  life-testing  experiments 
when  the  observations  do  arise  in  a  certain  order  and  it  is  possible  for  the 
experimenter  to  select  a  few  quantiles,  the  choice  of  the  number  of  quantiles 
and  the  spacings  between  the  quantiles  being  left  to  the  experimenter.  Ex¬ 
pression  for  the  quantile  estimators  of  the  parameters  are  discussed  in  this 
section.  The  scale  parameter  is  assumed  to  be  kner/n,  the  optimum  spacing 
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of  the  quantiles  for  estimating  the  locatio  i  parameter  y,  has  been  obtained 
for  any  number  of  quantiles.  If  the  location  parameter  (i  is  assumed  to  be 


known,  the  optimum  spacing  of  the  quantiles  foi  estimating  q  has  been  deter¬ 
mined  for  the  number  of  quantiles  =*  3,4. 

Let  x*  %  <  x 
v  - 

of  size  n  from  the  logistic  distribution  (1.1).  The  following  expressions 
will  be  needed  and  are  defined  as  follows 


(n2)^-^X(nk) 


be  the  k  order  statistics  m  a  sample 


(2.1) 

(2.2) 

(2.3) 

(2.4) 

(2.5) 


(2.6) 


Xi  *  n  »  Xo  -  °»  Xk+1  -  i  -  2,  ...,  I: 

n  -  » 

pUi 

\  -  J  f(t)dt,  UjL  »  log^/U-xJ),  i  »  1,  2,  ...,  k 


fi  *  f^ui^  “  xi(1_xi)»  fo  2  fk+i  i  *1.  2,  ...,  k 


k+1 


X*  X  ^i^i-l^M^i^n )  'H-l^^i-l^n.  )* 

i-1  1 


k+1 

Y*  I  C(xi(i-Xi)ioge(\i/(1-\i)) 
i«l 

■Xi.i(1-Xi.i)l°ge(Xi.i/(1-Xi.i))(Xi(1-Xi)X(ni) 

’Xi-l^1"Xi-l^x(ni_1)^+^xi"xi-l^ 


k+1 


*1  *  Z  ^xi"xi-l^x"xi"xi-i^ 


i»l 

^  (X1(l-Xt)loge(Xi/(l>Xi))  -  Xi.1(l-X1.1)logp(l..1/(l-\i_1))r 

^  Xi  ”  Xi-i 

i«l  1  11 


(2.7)  Kg 
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k+1 

(2.8)  Kg-  £  (l-xi-Xi.1)(\i^-Xi)loge(Xi/(l-\i)) 
i=l 

"  Xi-l(l‘Xi-l)l0MXi-l/(l'Xi-l^  • 


Three  different  cases  will  be  discussed. 

(A)  Estimation  of  u,  (q-  known) 

From  the  general  expressions  for  the  estimators  derived  by  Ogawa 
(1951) »  the  best  linear  unbiased  estimator  p*  of  p,  and  its  variance  are 
given  by 

(2.9)  «  jq  - 

(2.10)  V(„*)=2i  .  ^ 

(2.9)  and  (2.10)  give  the  estimator  of  p*  and  its  variance  for  a  fixed 
number  k  of  quantiles  and  for  fixed  values  of  For  a  fixed  number  of 

quantiles,  the  following  theorem  gives  the  values  of  X^ 's  which  will  minim 
ize  V(p*). 

Theorem  1.  For  a  fixed  number  k  of  sample  quantiles,  the  spacing  of  the 
quantiles  for  which  V(p*)  is  minimal,  is  symmetric  and  is  given  by 
\t  »  i/(k+l)  . 

Since  minimizing  V(u*)  would  be  equivalent  to  maximizing  then 

by  maximizing  K.  w.r.t.  X^s  (i  »  1,  2,  ...,  k),  the  optimum  Xj's  will 
be  obtained  as  the  solutions  cf  the  system  of  equations 


(2.11) 


Xi+1  '  Xi  "  Xi  "  Xi-1 


i  =  1,  2,  ..<) 


k 
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The  above  se^  of  equations  is  satisfied  for  Xi  =  i/(k+i).To  prove  that 

2 

this  set  of  X^’s  does  maximize  K^,  consider  the  matrix  D  •  (3  K^/dX^dX^) 
evaluated  at  X^  »  i/(k+l).Let  be  the  determinant  of  the  matrix  ob¬ 

tained  from  D  by  deleting  the  last  p  rows  and  p  columns.  Then 


■  -  2(Vi!  -  V>  *  • 

It  follows  that  the  matrix  D  is  negative  definite.  Hence  is 
maximized  when  «  i/(k+l).  Now  for  \±  *  i/(k+l)  and  for  i  =  *,  2,  ...,  k, 
it  follows 


(2.12) 


Xi  +  Xk-i+l  “  1 


ui  +  w  ■  ° 


fi  "  fk-i+l  • 


Clearly,  (2.12)  implies  that  the  spacing  is  symmetric  and  =  0. 
Thus  the  best  linear  unbiased  estimator  with  optimum  spacings  is  given  by 


(2.13) 


(2.14) 


“*  -  EBT-afotei  I  . 

1  * 


i=l 

v&i»)  - 

a<:nk(k+2) 


The  Craoer-Rao  lower  bound  for  the  variance  of  an  unbiased  estimator  for  t* 
and  cr  is  known  is  giver  by  3o2/(a2n). Hence  the  relative  efficiency  is 
given  by 
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^  C.R.  lower  bound  k(k+2) 

(  •  4  ^5  =  ^7  ■ 

The  above  relati  e  efficiency  increases  with  k,  its  minimum  (for  k  >  0) 
being  0.75. 

Estimation  of  u  for  censored  samples 

In  practice  there  often  arise  situations  when  sane  of  the  observations 
are  missing.  Suppose  that  the  (r^l)  smallest  and  the  (r2-l)  largest 
observations  are  not  available.  Then  imposing  the  following  restriction  on  k 

(2.16)  k  <  min  (^f-  , 

the  best  linear  unbiased  estimator  for  p  can  be  obtained  from  (2.13). 

(B)  Estimation  of  g  when  y,  is  known 

In  this  case  it  has  been  proved  by  Ogawa  that  the  best  linear 
unbiased  estimator  for  g,  for  a  fixed  number  k  of  sample  quantiles  is 

(2.17)  g*  =  Y/Kg  -  y.  Kj/Kg 

(2.18)  V(o*)  *  a2/(nKg) 

Kg  and  being  defined  by  (2.7),  (2.8), 

As  in  the  previous  case,  the  problem  of  deriving  the  optimum  spacing  of 


.  * 


the  sample  quantiles  arises.  The  optimum  spacing  for  g*  is  obtained  by 
maximizing  Kg.  Now 
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dK2 

d\.  "  °»  for  1  =  1>  2»  “•»  k 


yields 


(2.19) 


c 


Q, 


«i+1 


Xi"Xi-l 


Xi+l”Xi 


Q. 


i+i 


Xi'Xi-l 


Xi+l‘Xi 


"|  =  0,  (i  -  1,  2,  ...,  k) 


where 


(2.20)  2,  =  1,0  -h  >10*,  (  A-}  -  ^  ) 


For  a  subclass  of  the  class  of  all  distribution  function,  Tischendorf 
(1955)  has  derived  necessary  conditions  for  the  spacings  that  makes  Kg 
maximum.  It  can  be  verified  that  the  logistic  distribution  belongs  to  this 
subclass.  This  necessary  condition  for  logistic  distribution  is 


(2.21)  2{(i-2Xi)loge  (5^)  ,  l}  -  jf-A-  -  Atij-  0,  1,1,2,  ....  *  . 


From  (2.20)  and  (2.21)  it  is  clear  that  the  optimum  spacing  can  be  ob¬ 
tained  by  solving  (2.21).  However,  it  is  not  possible  to  solve  (2.21)  explic¬ 
itly  for  X Besides  the  system  of  equations  (2.21)  may  also  possess  multiple 
roots  which  further  entails  a  choice  of  the  proper  X^’s.  A  slight  simplifi¬ 
cation  of  the  problem  is  effected  by  considering  only  rymmetric  spacings. 

For  k  ■  2,  and  using  symmetric  quantiles, 


L 
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(2.22) 


2Xi(1-Xi)2(10(3ef\i/(1-\i)}  )2 
*2  - - - 


For  this  case,  equation  (2.21)  becomes 


(2.23) 


1  -  3XX  + 


J/Vi  ^  ^ 

l-2x1  loee  C  1-XX  > 


0  . 


By  solving  (2.23) >  it  was  found  that  =  *1°3  is  a  solution  of  (2.23)  and 
it  vas  verified  that  Kg  does  have  an  absolute  maximum  at  X^  -  .103  • 

The  estimator  a*  and  its  variance  are  given  by 


(2.24) 

(2.24) 


a*  =  .4l92(.\'£|.  #897n]+1)"x([.  3.03113+!)) 


V(o*)  =  1.0227  ~  . 


It  can  be  shor.m  that  the  Craner-Rao  lower  bound  for  the  variance  of  an 
unbiased  estimator  8  of  0  is  , 


(2.26) 


V( 6)  >  — S2-* 
n(3+rr) 


[The  details  are  lengthy  and  have  been  omitted.]  Hence  the  relative  efficiency 
of  a*  as  compared  with  the  Cramer -Rao  bound  is  68.38°/0  . 

For  k  »  3j  the  assumption  that  the  spacing  of  the  quantiles  is  symmetric 
gives  Xg  to  be  equal  to  .5>  and  the  coefficient  of  n]+l)  *n  the 

estimator  of  a  is  zero.  In  general  for  k  =  2ra  +  1,  the  condition  of 
symmetric  spacing  reduces  the  coefficient  of  n]+l)  *n  es*,i" 

nator  of  a  to  zero.  Thus  for  k  *  3,  Kg  is  the  same  function  of  X-^  as 


given  in  (2.22)  so  that  the  estinv  tor  a*  and  its  variance  is  again  given 
by  (2.24)  and  (2.25),  respective! i. 
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(C )  Estimation  of  both  ^  and  q 

In  the  case  where  both  ^  and  a  are  unkncr.m,  the  estimators  are 

given  by 


(2.27) 

•**  -  i  oyt  - 

(2.28) 

a*  =  -  (-K3X  +  KjY) 

where 

A  *  KjK^  -  . 

(2.29) 

Var(M,*)  =  2?  h.  ,  v(o*) 

2  Ko 

Cov  (n*,  or*)  »  -  ~  ~ 

The  problem  of  obtaining  the  optimum  spacing  in  this  case  even  under  the 
simplifying  assumption  of  symmetry  of  the  spacings  is  complicated  since  mini¬ 
mizing  the  generalized  variance  leads  to  simultaneous  equations  which  cannot 
be  solved  explicitly. 

3.  Blom's  and  Jung’s  Estimators 

This  section  discusses  Blom's  and  Jung's  estimators  for  estimating 
tt  and  a  by  linear  functions  of  order  statistics.  Blom  (1957)  approximated 
the  best  linear  unbiased  estimators  by  estimators  that  are  unbiased  but  do  not 


necessarily  have  the  minimum  variance  of  all  linear  unbiased  estimators. 

Jung  (1955)  approximited  the  best  linear  unbiased  estimators  by  estimators 
that  are  "asymptotically  unbiased  and  asymptotically  best”. 

Since  this  investigation  was  carried  out,  the  best  linear  unbiased 
estimators  of  the  location  and  scale  parameters  using  order  statistics  have 
been  computed  for  sample  size  <  25,  and  are  given  in  Gupta,  Qureishi  and 
Shah  (1965)*  However  the  estimators  of  Blom  and  Jung  are  relatively  simple 
to  compute  and  hence  it  is  of  interest  to  determine  how  good  these  asymptot¬ 
ically  'good'  estimators  are  when  compared  to  the  best  linear  unbiased  esti¬ 
mators  in  terms  of  their  relative  efficiency  for  moderate  sample  sizes. 

This  investigation  would  thus  be  useful  for  estimating  the  parameters  for 
sample  sizes  >  25,  and  also  as  providing  some  indication  of  the  general 
properties  of  these  two  kinds  of  estimators,  and  hence  has  been  included. 

The  authors  computed  both  these  estimators  for  both  \i  and  p.  Since 
the  estimators  of  Jung  are  biased  they  were  modified  by  multiplying  them  by 
an  appropriate  constant.  The  foilwing  were  the  main  results. 

(i)  The  modified  estimator  of  Jung  for  estimating  u  reduces 
to  Blom's  estimator  n*  where 

*  I  n(n'lJrnT2j  *  X(i) 

i«l 


The  approximate  variance  of  n1  is  given  by 


<  d 


Jung's  estimation  of  o  modified  to  make  it  unbiased  did  not 
reduce  to  Blon's  estimator  and  its  exact  variance  for  n  <  25 
was  found  to  be  uniformly  lower  than  the  exact  variance  of  the 
corresponding  estimator  of  Blom.  Jung's  estimator  for  c,  a 
is  given  by 


- v— — p -  •  Y  [-(n+l)2+  2«j(n+l)  +  2v(n+l-v)losp(v/n+l-v)]  x,  »  . 

n(n+l)  i3+rr  )/3  ^1 


Blom's  estimator  of  <r>  o'  is 


*  Ivin 


where 


ai(n+l-i)(ci-ci-1) 

d(n+l)2 


>  i  <  n 


i(n+l-i 

(n+1)2 


2 1 


[i+l) (n-i ] 
(n+l)2 


^(i+l.n),  i  <  n 


i  £  t  =  -  ^(n-i+l.n),  n  -  i  >  i  -  1 


l<  ■ 


Table  I  gives  the  coefficients  of  the  order  statistics  in  Blom's 
estimator  of  |*.  The  last  column  gives  the  exact  variance  of  the  estimator. 
These  exact  variances  were  calculated  by  using  the  variances  and  covariances 
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or  -the  order  statistics  given  in  Shah  (19&5)  and  Gupta,  Shah  and  Qureishi 
^1965)-  The  relative  efficiency  of  this  estimator  uhen  compared  with  the 
best  linear  unbiased  estimator  is  given  in  Table  III  for  selected  values  of  n. 

Table  II  gives  the  coefficients  of  the  order  statistics  in  the  modified 
Jung's  estimator  of  a,  since  this  estimator  has  smaller  variance  than  Blom's 
estimator.  Again  the  last  column  gives  the  exact  variances  of  these  estimators. 

Although  Blom's  estimator  of  o  did  have  higher  variance  than  the  corres¬ 
ponding  estimator  of  Jung  (modified),  yet  it  does  have  fairly  high  relative 
efficiency  when  compared  with  the  best  linear  unbiased  estimator,  and  was 
found  to  be  at  least  96  percent  for  n  <  25.  Table  III  gives  the  relative 
efficiency  of  Blom's  estimator  of  o,  and  the  relative  efficiency  of  Jung's 
estimator  (modified)  of  a  for  selected  values  of  n. 

From  this  table  it  is  obvious  that  Blom’s  and  Jung's  estimators,  besides 
being  simple  to  compute,  have  very  high  relative  efficiency  even  for  moderate 
values  of  n,  so  that  for  n  >  25,  one  could  expect  these  estimators  to  be 
almost  as  efficient  as  the  best  linear  unbiased  estimators  of  Lloyd. 
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Table  I 

Coefficients  of  the  ith  order  statistic  in  the  unbiased  nearly  best  estimator  cf  y 

tie  mean  of  the  logistic  distribution,  using  Blom's  method. 

Coefficient  of  x,  .  , .  =  coefficient  of  x,... 

(n-i+1)  U) 


\i 

n\ 

1 

2 

3 

4  5 

6789 

10 

11  12  13 

a 

5 

.1429  .2286 

.2571 

.1927 

6 

.1071 

.1786 

.2143 

.1594 

7 

.0833 

.1429  .1786 

.1905 

.1358 

8 

.0667 

.1167 

.1500 

.1667 

.1182 

9 

.0545  .0970 

.1273 

.1455  .1515 

.1047 

10 

.0455 

.0818 

.1091 

.1273  .1364 

.^939 

11 

.0385 

.0699  .0944 

.1119  .1224 

.1259 

.0852 

12 

.0330 

.c6o4 

.0324 

.0989  .1099 

.1154 

.0779 

13 

.0286 

.0527 

.0725 

.0879  .0989 

.1055  .1077 

.0717 

14 

.0250 

.0464 

.0643 

.0786  .0893 

.0964  .1000 

.0665 

15 

.0221 

.0412 

.0574  .0706  .0809 

.0882  .C926  .0941 

.0620 

16 

.0196 

.0368 

.0515 

.0637  .0735 

.0809  .0858  .0882 

.0580 

17 

.0175 

.0330 

.0464 

.0578  .0671 

.0743  .0795  .0826  .0836 

.0546 

18 

.0158 

.0298 

.0421 

.0526  .06l4 

.0684  .0737  .0772  .0789 

.0515 

19 

.0143 

.0271 

.0383 

.0481  .0564 

.0632  .0684  .0722  .0744 

.0752 

.0487 

20 

.0130 

.0247 

.0351 

.0442  .0519 

.0584  .C636  .0675  .0701 

.0714 

.0463 

21 

.0119 

.0226 

.0322 

.0407  .0480 

.0542  .0593  .0632  ,c66l 

.0678 

.0683  .0440 

22 

.0109 

.0208 

.0296 

.0375  .0445 

.0504  .0553  .0593  .0623  .0642 

.0652  .0420 

23 

.0100 

.0191 

.0274 

.0348  .0413 

.0470  .0517  .0557  .0587 

.0609 

.0622  .0626  .0401 

24  , 

.0002 

.0177 

.0254 

.0323  .0385 

.0438  .0435  .0523  .0554 

.0577 

.0592  .0600  .0385 

25  . 

.0085  .0164 

.0236 

.0301  .0359 

.0410  .0455  .0492  .0523 

.0547 

.0564  .0574  .0573  .0369 

Table  II 


Coefficients  of  the  (n-i+l)th  order  statistic  x(n_i+1'  *n  the 


linear  estimator  of 

a  (by  Jung's  method)  modified  to  make  it 

unbi 

ased. 

1  2  3 

4 

5 

6  7  8  9  10  11 

12 

13 

Variance 

2 

. 0 

5 

.3538  .2038  0 

.1706 

6 

.2907  .2024  .0715 

.1372 

7 

.2459  .1907  .1024 

0 

.1147 

8 

.2125  .1767  .U47 

.0396 

.0985 

9 

.1867  .1630.  .1180 

.0616 

0 

.0864 

10 

.1663  .1503  .1170 

.0737 

.0251 

.0769 

11 

.1497  .1389  .1138 

.0800 

.0412 

0 

.0693 

12 

.1360  .1288  .1095 

.0823 

.0514  .0174 

.0630 

13 

.1244  .1198  .1049 

.0834 

.0577 

.0295  0 

.0578 

l4 

.1147  .1119  .1001 

.0827 

.0615 

.0379  -0128 

.0534 

15 

.1062  .1048  .0955 

.0813 

.0636 

.0436  .0222  0 

.0496 

16 

.1989  .0934  .0911 

.0793 

.0645 

.0475  .0291  .0098 

.0463 

17 

.0925  .0927  .0869  -0771 

.0646 

.0500  .0341  .0173  0 

.0434 

18 

.0869  .0876  .0629 

.0748 

.0641 

.0516  .0377  .0230  .0077 

.0409 

19  .0818  .0829  .0793 

.0724 

.0633 

.0524  .o4o4  .0274  .0138  0 

.0386 

20  , 

.0774  .0737  .0753 

.0700 

.0622 

.0523  .0422  .0307  .0187  .0063 

.0366 

21  , 

,0733  .0749  .0726 

.0677 

.0609 

.0527  .0434  .0332  .0225  .0113  0 

.0348 

22 

.0697  .0714  .0696 

.0655 

.0596 

.0524  .0441  .0351  .0255  .0154  .0052 

.0332 

23 

.0663  .0682  .0668 

.0633 

.0582 

.0518  .0445  .0365  .0278  .0188  .0095 

n 

.0317 

24 

.0633  .0652  ,0642 

.0612 

.0568 

.0511  .0446  .0374  .0296  .0215  .0130  . 

,0043 

.0303 

r-5  . 

,0605  .0625  .0618 

.0592 

.0553 

.0504  .0445  .0381  .0310  .0236  .0159  . 

0080 

n 

.2903 

-Computed  by  using  the  same  approximate  covariance  matrix  a.  used  in  Blom's  method. 


Table  III 


Table  giving  the  relative  efficiency  of  Blom's  estimator  of  and  Jung's 
estimator  (modified)  of  oj  relative  efficiency  is  with  refpect  to  the  best 
linear  unbiased  estimators. 


P.el. 

5 

7 

10 

15 

20 

25 

Blom's  Estimator 

of  li 

•991 

.993 

.996 

.997 

.998 

.999 

Jung's  Estimator 
(modified)  of  a 

•993 

•  993 

.999 

1.000 

1.000 

1.000 
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